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ANOTHER  POWER  IDENTITY  INVOLVING  BINOMIAL  THEOREM  AND 

FAULHABER’S  FORMULA 

PETRO  KOLOSOV 


Abstract.  In  this  paper,  we  derive  and  prove,  by  means  of  Binomial  theorem  and  Faulhaber’s 
formula,  the  following  identity  between  m-order  polynomials  in  T 

t  m  ra 

J2J2A^j(T-k)j  =  ^(-1  )m~kUm{i,k)-Tk  =  T2m+1,  £  =  T  <=  N. 

k= 1 j= 0  k= 0 


Contents 

1.  Introduction  and  Main  results 
1.1.  Error  of  approximation 

2.  Acknowledgements 

3.  Conclusion 

4.  Supplementary  Files 
References 


J 

13 

19 

19 

19 

19 


1.  Introduction  and  Main  results 
The  following  identity  holds  in  m-order  polynomials  Pm(£,T) 

i  m  m 

(1.1)  Pmi^T)  :  =  -  k)j  =  k)  ■  Tk  =  T2m+\ 

k= 1  j= 0  k= 0  e=Te  N 

where  £  G  Ni,m  G  Ni  and  T  G  M  and  Um(£,k),  Amj  are  coefficients  dependent  on  £  and  m, 
respectively.  The  term  (— l)m~^  appears  in  r.h.s  of  the  expression  (11.11)  by  means  of  sign-changing 
inside  J2jLo  Am.j  (T  —  k)i  by  Binomial  expansion  of  (T  —  k)P 

We  start  our  discussion  concerning  the  identity  (11. ID  from  the  derivation  of  the  partial  case  of 
the  polynomial  Pm(£,T)  for  m  =  1,  by  means  of  Faulhaber’s  identity  for  n3.  The  terms  of  the  of 
the  polynomial  P\  (T.  T)  =  )Cfc=i  A\:Jk3  (T  —  k'y  =  T3  over  k  produce  the  T-th  row  of  the 

sequence  A287326,  starting  from  k  =  1,  [3].  To  derive  the  partial  case  of  the  polynomial  Pm{£.  T) 
for  m  =  1,  let’s  consider  the  Faulhaber’s  identities  |TJ  for  odd  powers  n2m+1,  m  G  No 

=  0 

=  an1)  +  (?) 

=  120("+2)  +30(“J1)  +  (“) 

=  E  (2fc-l)!T(2n!,2fe)("+U1) 

1  <k<m 
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The  coefficients  T(m,  k )  in  these  formula  are  related  to  what  Riordan  [T9]  has  called  central  factorial 
numbers  of  the  second  kind.  In  his  notation, 

xm  =  T(m,k)x =  x(x  +  |  —  l)(x  +  |  —  2)  •  ■  ■  (x  +  |  +  1) 

1  <k<m 

The  coefficients  T(2m,2k)  are  always  integers,  because  the  adfc+2l  =  x^(x2/k4)  implies  recurrence 

T(2m  +  2,  2k)  =  k2T(2m,  2k)  +  T(2m,  2k  -  2). 

Central  factorial  numbers  of  the  second  kind  can  be  calculated  as  follows: 


(1.3) 


(2k  -  l)!T(2n,  2k)  =  -  £(- 1)>  ^ )  (r  -  j)2n , 
r  j= o  V  3  2 


where  r  =  n  —  k  +  1 .  The  formula  (11.31)  was  derived  by  Peter  Luschny  in 

The  Faulhaber’s  formula  and  polynomial  CCD  could  be  connected  the  following  way 

m+ 1 


n2m+l  =  J^(2k~  1)!T(2m  +  2,  2k)  lU  +  k  1 


k= 1 
n  m 


2k -l 


n 


n 


=  E  E  A™,jkj(n  -  ky  =  E;(-l )m-kUm(n,  k) 

k= 1 j=0  fc=0 

n2m~l  =  jT(2k  ~  l)\T(2m,2k)(n  +  ^  ^ 

k= 1  ' 

n  m— 1  m—  1 

=  E  E  Am-i^(n  -  ky  =  ^2  (~^)rn  l~kUrn-i(n,  k) 

k= 1  j= 0  k= 0 

The  forward  finite  difference  of  odd  power  could  be  reached  introducing  the  r  €  Ni  to  the  lower 
index  of  binomial  coefficient  of  (11.21)  the  following  way 

(r  =  1  :  An1  =  (©) 

=  6  ©;)  + 

A«5  =  120©,2)  +  30(5+;)  +  (©) 


r  =  1 
r  =  1 


=  0 

=  6("k)  +  (S) 

=  120(”J2)  +  30("JX)  +  (5) 


r  =  l:  An“-‘=  E  (2fc  -  l)!T(2m,  2*)(”+‘1l‘r)  =  E  (2*  -  l)!T(2m,2i)  (“+*-© 

1  <k<m 


-,2m— 1  _  fOU  OU\(  n+k— 1 

1  <k<m  _ 

I ! v  the  dentity  A f(x  —  1)  =  V/(x),  h  =  1,  backward  differences  could  be  reached  as  well 

fr=l:  W  =  (©)  =(V) 

V-3  =  6(3"r)  +  GZl)  =  6©  +  (y1) 

Vn5  =  120©;)  +  30©,,)  +  ©))  =  120("+‘)  +  30©  +  (©) 


r  =  1 
r  =  1 


r  =  1  :  Vn2”-^  £  (2fc  -  l)!T(2m,  2k)  (2"+Er)  =  E  (2*  -  l)!T(2m,  2/c)  (n+_2  ) 

l<fc<m  l</c<m 

Continuing  similarly,  for  every  r  >  1  we  get  Arn2m_1  and  Vrn2m_1.  To  derive  the  partial  case 
of  polynomial  Prn (£  =  T,T)  =  T2m+1,  T  €  Ni,  for  m  =  1,  recall  the  Faulhaber’s  identities 
An3  =  6 2 2 ' )  +  (o)>  Vn3  =  6(2)  +  ("o1).  thus,  the  perfect  cube  in  T  is: 

r3=x;iAT3(fc)  =  x;i6(fc+1)  +  w 

(1.4)  fc=0  fc=0  V  “ 

=Evm  =  E6f^  +  A'“1N 


k= 1 


k= 1 


0 
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Let’s  rewrite  the  expression  (II. 4|).  taking  to  the  attention  the  following  identity  (n^‘1)  =  0  +  1  +  2  + 
■  ■  ■  +  n,  thus 

(1.5)  T3  =  (1  +  6  •  0)  +  (1  +  6  •  0  +  6  •  1)  +  •  •  •  +  (1  +  6  •  0  +  •  •  •  +  6  •  (T  -  1)) 

Factorising  the  expression  (11.51)  we  get 

(1.6)  T3  =  T  +  (T  -  0)  •  6  ■  0  +  (T  -  1)  •  6  •  1  +  ■  ■  ■  +  (T  -  (T  -  1))  •  6  •  (T  -  1) 

Let’s  apply  the  sigma  notation  on  the  expression  (1 1 . 6 .  we  have 

T— 1  T— 1 

T3  =  T  +  ^  6  k{T -k)  =  J2  6fc(T  -  k)  +  1 

^  ^  k= 0  k= 0 

=  t + y  6fc(T  -k)  =  J2  6k(T  - k)  + 1  =  pi(T>  r)> 

k= 1  k= 1 

see  (11.11)  for  P\(T,T). 

Although,  the  expression  (11.71)  is  derived,  essentially,  from  identity  in  finite  differences  of  cubes, 
such  as  T3  =  YlkZ, o  A T3(k)  and  T3  =  Y^=i  VT3(fe),  the  finite  differences  AT3(/c)  /  6k(T  —  k)  +  1 
and  VT3(fc)  /  6 k(T  -  k)  +  1. 

Over  the  manuscript  we  review  the  Pm(£,T)  in  sense  of  backward  difference  VT2m+1,  m  £  Ni, 
i.e  we  use  the  summation  limits  over  k  as  1  <  k  <  £  in  dm 

Therefore,  we  have  reached  the  generating  function  of  A287326  in  r.h.s  of  (jl.7|> .  The  corre¬ 
sponding  coefficients  Amj  in  the  dehnition  Aijk^  (T  —  k)i  =  1  +  6 k(T  —  k )  of  generating 

function  are:  Ago  =  1,  Agi  =  6.  Note  that  by  definition  the  polynomials  Am ,j k3  {T  —  k)J 

should  be  displayed  starting  from  the  Amfik°(T  —  k)°,  we  will  denote  these  polynomials  starting 
from  the  m-th  term,  as  the  order  of  summation  doesn’t  change  it’s  result.  Let’s  construct  the 
triangle  A287326,  every  (£  =  T)-th  row  sum  starting  from  k  =  1  gives  the  terms  of  polynomial 
Pi(£,T)  =  Ylk=i  X^=i o-A-ijk^T  —  k)i  =  J2k=i  §k(T  —  k)  +  1,  starting  from  k  =  1.  We  review  the 
triangle  A287326  in  context  of  the  partial  case  of  polynomial  Pm(£,T),  m  =  1,  £  =  T  G  No 


£  =  T  =  0 
£  =  T  =  1 

£=T= 2  1 

£  =  T  =  3  1 

£  =  T  =  4  1  19 

£  =  T  =  5  1  25 


1 

1  1 

7  1 

13  13  1 

25  19  1 

37  37  25  1 


Table  1.  Triangle  generated  by  the  polynomial  Xg=o  AigAd  (T  —  fc)J  =  6 k(T  —  k)  +  1,  0  < 
k  <  t  =  T  S  No,  sequence  A287326  in  OEIS.  Summation  of  the  {£  =  T)-th  row  terms  from 
k  =  1  gives  P1(T,T)  =  T3,  see  (11.11). 


The  sum  of  the  (£  =  T)-th  row  terms  of  Itablc  T1  starting  from  k  =  1  generates  the  partial  case 
of  (11.11)  for  to  =  1  and  £  =  T  €  N,  that  is  Pi(T,  T)  =  T3.  Binomial  distribution  of  the  row 
terms  ofltable  Tlcan  be  easily  proven  by  reviewing  of  its  generating  function  A\jkP  (T  —  k)i  = 

6k(T  —  k)  +  1  =  6kT  —  6 k2  +  1,  which  is  parabolic  for  every  given  variable  T,  and,  therefore,  is 
symmetrical  over 

Let’s  show  a  few  properties  of  the  generating  function  of  P\(T,T),  define  the  term  D\(T,k ) 
(1.8)  L>i(T,  k)  :=  Ahikj(T  -  k)j  =  6k(T  ~k)  +  l 

l<j<0 

The  following  properties  hold  in  terms  of  D\(T,  k ) 
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(1)  Finite  difference  property,  A (n3)  =  D\  ^nZ+2n+2, 1^ ,  where  n2+™+2  is  the  n-th  Central 
polygonal  number,  see  [22]. 

(2)  Binomial  transform  of  the  second  order  of  D\(n,  k):  Di(n  +  l,k)  =  2D\(n,  k)  —  Di(n—  1,  k). 

(3)  Symmetry,  Di(n,k)  =  D\{n,n  —  k),  0  <  k  <  n. 

Below  we  show  initial  ten  polynomials  Pm(£,T)  =  'Y^k=o{^)m~k^rn{^k)  ■  Tk,  generated  by  the 
partial  case  of  (11.111  for  m  =  1  over  £,  0  <  £  <  10: 


(1.9) 


l  1  L  1 


E  6 k(T  -  k)  +  1 

=  £(-+ 

-fec/i(4fc)  -Tk 

k=  1  j= 

=0 

fc=i 

k= 0 

'£  = 

1 

-5  +  6  T 

=  -C/i(l,0)  - 

'j'O 

+  I7i(l,l) 

T 1  = 

£  = 

2 

-28  +  18T 

=  — Ci(2,  0) 

'j'O 

+  U\(2, 1) 

T1  = 

-Pi  (2,  T) 

£  = 

3 

-81  +  36T 

=  — C7i  (3, 0) 

'j'O 

+  Ui(3, 1) 

T1  = 

-£4(3,  T) 

£  = 

4 

-176  +  60T 

=  -C7i(4,0)- 

'j'O 

+  £4(4,1) 

T1  = 

-£4(4,  T) 

£  = 

5 

-325  +  90T 

=  — C/i(5, 0) 

'j'O 

+  £/i(5,l) 

T1  = 

Pi&T) 

£  = 

6 

-540  +  126T 

=  -C/i(6,0)  - 

'j'O 

+  E£i(6,l) 

T1  = 

-£4(6,  T) 

£  = 

7 

-833  +  168T 

o' 

+ 

tS1 

1 

II 

'j'O 

+  C7i(7, 1) 

T1  = 

-£4(7,  T) 

£  = 

8 

-1216  +  216T 

o' 

OO 

tS1 

1 

II 

rpO 

+  £4(8,1) 

T1  = 

Pi(8,T) 

£  = 

9 

-1701  +  270T 

=  —Ui(9,  0) 

'j'O 

+  C+i  (9, 1) 

T1  = 

£4(9,  T) 

J  = 

10 

:  -2300  +  330T 

=  — C7i(10, 0) 

■T 

0  +  £4(10, 1)  • T 1 

=  £4(10 

T3  as 

T 

-+L 

The  coefficients  U\(£,k),  0  <  k  <  1  in  (11.91)  are  terms  of  the  sequence  A320047,  [21],  One  would 
ask,  ’’Why  to  show  the  identity  (11.11)  in  terms  of  (£  =  T)-th  row  ofltable  II  we  used  the  condition 
£  =  T  £  No,  but  to  show  the  identity  (11.11)  in  terms  of  (11.91)  we  apply  the  condition  T  — >  £,  and 
£  £  N?”  -  We  used  the  condition  T  — >  £,  and  £  £  N  to  show  the  identity  (11.11)  in  terms  of  (11.91)  as,  by 
definition  T  £  M  in  (11.11)  and  the  variable  £  £  N,  since  the  coefficient  Um(£,  k )  is  dependent  on  £,  and, 
therefore,  £  is  always  Ni,  so  if  we  speak  about  approximation  of  monomial  T2m+1  by  polynomial 
'  Tk  in  some  particular  point  £  £  N,  we  say  that  T  — >  £.  The  binomials,  listed 
in  (11.91)  are  linear  approximations  of  monomial  T3  in  some  neighborhood  of  the  point  £  £  N  and 
T  — >■  £.  The  following  figure  graphically  shows  the  binomials  J2k=o(~ ^-)1~kUi(£,  k)  ■  Tk ,  £=  1,2,3 
among  the  cubic  curve  in  T 
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Figure  1.  Linear  approximations  of  the  monomial  T3  by  binomials 
Pi(l,T),  Pi(2,T),  Pi(3,T)  in  some  neighborhood  of  the  point  l  =  1,2,3  with  T  — >  l,  see 

cm 

Now,  we  have  derived  and  discussed  the  pattern  A287326,  generated  by  the  E}=o  A\jk?(T  — 
ky  =  6k  (T  —  k)  +  1,  0  <  k  <  £  =  T  £  No-  Let’s  find  analogs  of  the  pattern  A287326  for 
m  >  1.  To  find  such  analogs,  let’s  determine  the  coefficients  A2tQ,  A2ji,  A22  ha  the  polynomial 
Ej=0  A2jkj(T  -  ky,  such  that  P2(i,  T)  =  ELi  Ej=0  A2jkj(T  -  ky  =  T5,  as  £  =  T  £  N,  by  the 
(1.2).  Firstly,  to  determine  the  coefficients  coefficients  A2iq,  A2,h  A2,2 ,  let’s  rewrite  the  polynomial 
P2  (£,  T )  =  Ek=i  Ej= 0  A2,jkj(T  —  ky ,  in  extended  view  as  follows 

EE'VAP’-fy 

fc= 1  j= 0 


(1.10) 


a2,2  k2(T-k)2  +  A2tl  k1(T-k)1  +  A2: 0  k°(T  —  k)° 

1 <k<£  1 <k<£  1 <k<£ 

A2,2  J2  k2(T2  -2Tk  +  k2)  +  A2:1  kT-k2  +  A2fi  Y  k\T  -  kf 

1 <k<e  1 <k<i  1 <k<e 

A2,2  Y  k2T2  -  2Tk3  +  k4  +  A2A  kT-k2  +  A2fi  Y  1 


l<k<£ 


l<k<£ 


l<k<£ 


aX2t2  y  k2  -  2A^T  E  +  A2,2  Y  k4  +  A^T  E  k 


,1  <k<£ 


,1  <k<£ 


.  l<k<£ 


,1 <k<£ 


-As,  1  E  k2  +  A2fi  E  1  =  T\  as  £  =  T  £  N. 


.  Kk<£ 


. l<k<£ 


Above  derivation  could  be  generalised  in  terms  of  Binomial  coefficients  and  Faulhaber’s  sum  for 
every  m  G  No  as  follows 


m  m—s 


£E(- lyAm.r. 

s= 0  j= 0 


m  —  s 


'  *  \ 

^  ^  Jtm+s  J  'J'm+s 


1 
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Note  that  the  part  A2io  Yli<k<i  k°(T  —  k)°  of  (11.1011  gives  an  indeterminate  form  as  k  =  T  since  the 
A2)o(T  —  T)°k°  contains  the  term  (T  —  T)°  =  0°.  Some  textbooks  leave  the  quantity  0°  undefined, 
because  the  functions  x°  and  0*  have  different  limiting  values  when  x  decreases  to  0.  For  our 
purposes  we  will  use  the  convention: 

Vx  :  x°  =  1, 

as  it  is  a  common  agreement,  see  [16].  Above  we  have  derived  an  expression  containing  sums  of 
powers  of  successive  natural  exponents,  where  the  powers  are  {1,2,3, 4}.  These  formulae  contain 
so-called  Bernoulli  numbers,  [7j.  By  Faulhaber’s  formula,  the  sums  of  successive  powers  in  k  are 
following 

e+i 


(Ml) 


£  * 

1  <k<£ 

£  *2 
1  <k<£ 

£  k3 

1  <k<£ 

£  k‘ 
l<k<e 


2£3  +  3£2  + £ 
6  ’ 

£4  +  2£3  +  £2 


6£5  +  15£4  +  10£3  -  £ 
30  ' 


Next,  we  substitute  the  identities  (11.1111.  into  (11. 1011.  respectively 


a2,2t 2 


(1.12) 


2£3  +  3£2  +  i 
6 


2A22T 


+A2^T 


l4  +  2£3  +  e 


+  A: 


2,2 


6£5  +  15£4  +  lOf5  -  £ 


30 


+  ' 


—  A2a 


2£3  +  3£2  +£ 
6 


+  A2fi£. 


Factorising  the  expression  (11.1211  and  collecting  the  terms  under  common  divisor  with  set  £  =  T  €  N, 
we  get 

_  rp 


(1.13) 


A2t2T5  —  A2,2T  +  30A2jo  -l  —  -l  5 
- 1-  2*2,1 - 7 —  =  1 


30  ’  6 

In  order  to  satisfy  (11.131)  for  each  £  =  T  £  N,  coefficients  A2^,  A2q,  A2)2  should  be  a  solutions  of 
following  system  of  equations 


(1.14) 


A 


2,1 


=  A 
=  1, 


30A2jq  —  A2  2  —  0. 


2,2 


The  solutions  to  the  system  (11.1411  are  following:  A2  2  =  30,  A21  =  0,  A2  0  =  1.  Hereby,  polynomial 
Y3= q  A2)jki(T  —  k)'J  takes  the  form 


(1.15)  A2,jkj(T  -  k)j  =  30k2  (T  -  k  f  +  1. 

3=0 

Let’s  construct  the  triangle,  every  (£  =  T)-th  row  sum  starting  from  k  =  1  gives  the  terms  of 
polynomial  P2{£,T )  =  Y%,=\  )Cj=o  A2.jkj(T  -  k)j  =  YX=i  30k2  (T  -  k)2  +  1,  over  k  in  range 
1  <  k  <  £  =  T  GNo-  We  review  the  triangle  in  context  of  the  partial  case  of  polynomial 
Pm(£,T),  m  =  2,  £  =  T€N0 
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0 

II 

bn 

II 

1 

£  =  T  =  1 

1 

1 

CN 

II 

bi 

II 

1 

31 

1 

£=T= 3  1 

121 

121 

1 

II 

II 

1 — 1 

271 

481 

271 

1 

£  =  T  =  5  1  481 

1081 

1081 

481  1 

Table  2.  Triangle  generated  by  the  polynomial  22 
k  <  £  =  T  £  No,  sequence  A300656  in  OEIS,  [8]. 
from  k  =  1  gives  P2(P, T)  =  T5,  see  (11.11). 

i=o  A2,jkj(T-k)j  =  30/c2(T— /c)2  +  l,  0  < 
Summation  of  the  (£  =  T)-th  row  terms 

By  the  identity  ED  in  polynomials  Pm(£,T),  the  sum  of  the  (T  =  £)-th  row  terms  of  the 
Itable  21  over  k  from  1  to  £  =  T  £  N  gives  P-2(T,T )  =  T5.  Below  we  show  initial  ten  polynomials 
Pm(£,T)  =  l)m_fet/TO(£,  k)  ■  Tk,  generated  by  the  partial  case  of  (11.11)  for  m  =  2  over 

£,  1  <  £  <  10: 

(1.16) 

e  2  t  2 

p2(i,  =  A^kj(T  -  ky  =  30fc2(T  - k)2  + 1  =  fc)  • 

fc=l  j=0  k=  1  fc=0 


'£  =  1  :  31-60T  +  30T2 
£  =  2  :  512  -  540T+  150T2 
£  =  3  :  2943  -  2160T  +  420P2 
£  =  4  :  10624  -  6000T  +  900T2 
£  =  5  :  29375  -  13500T  +  1650T2 
t  =  6  :  68256  -  26460T  +  2730 T2 
£  =  7  :  140287  -  47040T  +  4200T2 
£  =  8  :  263168  -  77760T  +  6120T2 
£  =  9  :  459999  -  121500T+  8550T2 
£  =  10  :  760000  -  181500T  +  11550P2 


=  U2(l,  0)  •  T°  -  U2{  1, 1)  •  T1  +  C/2 (1,  2)  •  T2  =  P2(l,  T) 

=  U2{ 2,  0)  ■  T°  -  £/2(2, 1)  •  T1  +  C/2 (2,  2)  •  T2  =  P2(2,  T) 

=  C/2 (3,  0)  ■  T°  -  C/2(3, 1)  •  T1  +  0/2(3, 2)  •  T2  =  P2(3,  T) 

=  0/2(4,  0)  •  T°  -  £/2(4, 1)  •  T1  +  U2{ 4,  2)  •  T2  =  P2(4,  T) 

=  0/2(5,  0)  •  T°  -  0/2(5, 1)  •  T1  +  0/2(5,  2)  •  T2  =  P2( 5,  T) 

=  C/2 (6, 0)  ■  T°  -  0(2(6, 1)  •  T1  +  U2{ 6, 2)  •  T2  =  P2(6,  T) 

=  C/2(7,  0)  •  T°  -  0/2(7, 1)  •  T1  +  0/2(7, 2)  •  T2  =  P2( 7,  T) 

=  0/2(8,  0)  •  T°  -  0/2(8, 1)  •  T1  +  0/2(8,  2)  •  T2  =  P2( 8,  T) 

=  0/2(9,  0)  •  T°  -  0/2(9, 1)  •  T1  +  0/2(9, 2)  •  T2  =  P2(9,  T) 

=  0/2(10, 0)  ■  T°  -  C/2(10, 1)  ■  T1  +  0/2(10, 2)  ■  T2  =  P2(10,  T) 


=  T5  as  T  -a  £. 


The  coefficients  0/2(0,  Oc),  0  <  A;  <  2  in  (11.161)  are  terms  of  the  sequence  A316349,  [I?].  Above  poly¬ 
nomials  are  approximations  of  monomial  T5  in  some  neighborhood  of  the  point  T  — >  £,  the  following 
graph  shows  corresponding  polynomials  22k=i  ^.?=o  A'2,jk] (T—ky  =  ^2=0(— l)2~fct/2(f(,  k)-Tk,  £  = 
1, 2, ..  and  the  T5  curve 


PETRO  KOLOSOV 


Fen 


Figure  2.  Approximations  of  monomial  T 5  by  polynomials  P2(1,T),  P2(2,T),  P2(3,T)  in 
some  neighborhood  of  the  point  ^  =  1,2, 3,  with  T  — >  i  see  (11.161). 


Similarly,  finding  the  coefficients  A3>0,  A3> i,  A3j2,  A3)3  in  P3(£,T )  =  Y?k=i  X^=o  A3,jkj(T  -  k)j , 
we  get  A3io  =  1,  A3ti  =  0,  A3< 2  =  — 14,  A3)3  =  140,  therefore,  for  each  £  =  T  €  N: 

T7  =  140/c3  (T  —  A;)3  —  14/c2(T  —  k)2  +  1. 

1  <k<i 

Below  we  show  a  few  initial  rows  of  the  triangle,  every  {£  =  T)-th  row  sum  starting  from  k  =  1  gives 
the  terms  of  polynomial  P3(£,  T)  =  Ylk=i  Sj=o  7l3j£:J(T  —  =  Ylk= l  140A;3(T  —  k)3  —  14k2 (T  — 

k )2  +  1,  over  k  in  range  1  <  k  <  £  =  T  (L  No-  We  review  the  triangle  in  context  of  the  partial  case 
of  polynomial  Pm(£,T),  m  =  3,  £  =  T  £  No 


£  =  T  =  0  1 

£  =  T  = 1  1  1 

£  =  T  =  2  1  127  1 

£  =  T  =  3  1  1093  1093  1 

£  =  T  =  4  1  3793  8905  3793  1 

£  =  T  =  5  1  8905  30157  30157  8905  1 

Table  3.  Triangle  generated  by  the  polynomial  Y^j=o  A3jki(T  —  k)i  =  140 fc3(T  —  k )3  — 
14fc2(T  —  k )2  +  1,  0  <  k  <  £  =  T  £  No,  sequence  A300785  in  OEIS,  [9j.  Summation  of  the 
(£  =  T)- th  row  terms  from  k  =  1  gives  P3(T,  T)  =  T7,  see  (11.11). 


As  identity  (11.11)  holds,  the  sum  of  the  (T  =  ^)-th  row  terms  of  Itable  ~3l  over  k  from  1  to  T  =  £ 
equals  to  T2m+1,  m  =  3.  For  the  case  P3(£,T),  the  following  polynomials  Yll=i  A3.3£A {T  — 
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k)i  =  Ylk=o(~ l)3  kU3(i,  k)  ■  Tk,  for  £  =  1,2, 10  can  be  generated 
(1.17) 

l  3  3 

P3(£,T)  =^2^2  a 3,jkj{T  -  k)j  =  J2  140fc3(T  -  k)3  -  Uk2(T  -  k )2  +  1  =  ^(-l)3-fcF3(£,  Jfe)  •  Tk 
fe= i j=o  i <k<e  fe= o 

'  1  =  1:  —125  +  406T  —  420T2  +  140T3  =  P3(1,T) 

£  =  2  :  -9028  +  13818T  -  7140T2  +  1260T3  =  P3(2,T) 

£  =  3  :  -110961  +  115836T  -  41160T2  +  5040T3  =  P3(3,  T) 

£  =  4  :  -684176  +  545860T  -  148680T2  +  14000T3  =  P3  (4,  T) 

£  =  5  :  —2871325  +  1858290T  —  411180T2  +  31500T3  =  P3(5,  T) 

"  '  £  =  6  :  -9402660  +  5124126T  -  955500 T2  +  61740T3  =  P3(6,  T) 

£  =  7  :  -25872833  +  12182968P  -  1963920 P2  +  109760T3  =  P3( 7,  T) 

£  =  8  :  -62572096  +  25945416T  -  3684240 T2  +  181440T3  =  P3  (8,  T) 

£  =  9  :  -136972701  +  50745870T  -  6439860T2  +  283500T3  =  P3(9,  T) 

£=ID:  -276971300  +  92745730T  -  10639860T2  +  423500T3  =  P3(10,  T) 

=  T7  as  T  £. 

The  coefficients  U3(T.  k),  0  <  k  <  3  in  (J1.17I)  are  terms  of  the  sequence  A316387,  [18].  Above 
polynomials  are  approximations  of  monomial  T‘  in  some  neighborhood  of  the  point  T  — >  t.  The 
following  graph  shows  corresponding  polynomials  Ylk=o(~  1)3— fcC/3 (P,  k)  ■  Tk,  T  =  1,2,3  and  the 
T‘  curve 


Figure  3.  Approximations  of  monomial  T7  by  polynomials  P3(l,  T),  P3(2,T),  P3(3,T)in 
some  neighborhood  of  the  point  £=1,2, 3,  with  T  — >  £,  see  (11.171). 

Now,  let  generalise  the  the  identity  Pm(T,T )  =  T2m+ 1 ,  T  €  N  for  every  integer  m  >  0,  as  we 
have  already  discussed  the  cases  m  =  1,2,3.  Firstly,  let’s  find  a  recurrence  formula  for  coefficients 
Amj,  0  <  j  <  m,  see  definition  (11.11).  To  reach  a  recurrent  formula  of  Amj,  m  >  0,  first  let’s  fix 
the  unused  values  of  Amj  =  0,  for  j  <  0  or  j  >  m,  so  we  don’t  need  to  consider  the  summation 
range  for  j.  Let’s  rewrite  the  polynomial  Pm(£,  T)  again 

t  m  m 

Pm(£,T)  =  =  ^2(-l)m-kUm(£,  k)  -Tk, 

k= 1  j= 0  k= 0 


(1.18) 
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Since  the  symmetry  holds 


Dm(T,  k )  :=  J2  Am,jkj(T  -  k)j ,  Dm(T ,  k)  =  Dm(T  —  k,  k) 

3=  o 

We  can  rewrite  (11.1811  as  follows 

£  m  t—  1  m 

(1.19)  E  E  Am,Jkj  ( T  -  k)j  =  E  E  A^kj  ( T  -  ky 

k= 1  j= 0  fc=0  j= 0 

By  expanding  (T  —  k)J  in  r.h.s  of  ()  1 . 1 9 1)  and  using  Faulhaber’s  formula  [IT],  the  result  is 
e-i  i-  i 

'3 

A. 


E<r  -  fc)3V  =  EE  ( i  pN-i)v^ 


k=0 


k= 0  i 


(1.20) 


.  \i)  i  +  j  +  1 


E 

_  t 


i  +  j  +  1 
t 


y-  A\  (-1)'  A  +  J  +  1 
^  w  *  +  J  + 1  V  t 


Btli+j+ 1-*  -  Si+i+i 


'.A  (-1)' 


BfTj-ifi+j+i-i  _  ^  I  J  _ 


V  *  +  j  +  1 


Si+j+iT^-* 


(*)  W 

where  Bt  are  Bernoulli  numbers  [IT].  Now,  we  notice  that 

"A  (~A  (i+3  +  l\=i  (2j+i1)(2/) ’  lf  1  =  0; 


(1.21) 


\ij  z  +  j  +  1  V  t 


In  particular,  the  last  sum  is  zero  for  0  <  t  <  j.  Now,  by  substituting  to  the  (*)  part  of  (11.201)  the 
result  of  (11.211).  we  have 


(1.22) 


i,t 


A  (-A  fi+j  +  i 


^  '  i)  i  +  j  +  1  V  i 


Rp-if+ih-t  _ 


(2j  +  l)(2/) 


TO  13+ 1 


+£ 


t>0 


(  A  |  j  )  BtT^~ili+^+ 1_t 

t  \2j -t  +  l) 


By  means  of  (11.221).  expression  (11.201)  takes  the  form 


£-1 


y  (t  -  kykj  = - t* ^-THj+i + y 

h  w+K1)  A 


(-It  (  3 


t  \2  j  —  t  +  l 


B  Tj-ifi+j+l-i 


(1.23) 


(*) 


Eh 


A  (-A 


i)  i  +  j  +  1 


Bi+j+1T 


(O) 


We  have  to  remember  that  if  the  sum  over  some  variable  i  contains  (1),  then  instead  of  limiting 
its  summation  range  to  i  G  [0,  j],  we  can  let  i  G  [— oo, +oo]  since  (I)  =  0  for  i  outside  the  range 
i  G  [0,  j]  (i.e.,  when  i  <  0  or  i  >  j).  It’s  much  easier  to  review  such  sum  as  sum  from  — oo  to 
+oo  (unless  specified  otherwise),  where  only  a  finite  number  of  terms  are  nonzero,  this  fact  is  also 
discussed  in  [12] .  Therefore,  we  haven’t  shown  detailed  bounds  of  summation  in  above  derivation. 
Now,  we  keep  our  attention  on  (|1 ,23|) .  To  combine  or  cancel  identical  terms  across  the  two  sums 
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in  (11.2311  more  easily,  let  introduce  k  =  2j  +  1  —  t  to  (*)  and  k  =  j  —  i  to  (o) ,  respectively 
i- 1 


J^(T  -  k)jkj 
=o 

£ 


1 


k= o  (2  J  +  1)(2/) 

(-i y 


THj+i  +  y — 
21 


(-IP  Pj 


(1.24) 

-V/  , 

kkJ  2j  +  1  —  k 

Let  be  £  =  T  in  (11.241).  thus 

T— 1 


2j  +  1  —  K  \K 


hTJ  'VK-'  j 


TD  rpK, 

£>2j+l-K-L 


J2(T~k)jkj 


k= 0 


(2j  +  i)(7) 


_r2j+l  V-  _ 
Z^  91 


(-IP  (j 


2j  +  1  —  K  \K 


TD  'T'K 

£>2j+l-K-L 


(1.25) 


E 


A  (-ip- 


B2j+i-KT* 


kJ  2  j  +  1  —  « 

= _ 1 _ 7"2i+1  _i_  2 

w +i)(?)  otr« 

Now,  using  the  definition  of  Amj, 

i- 1 


(-IP  (3 


2j  +  l  —  K\K 


TD  rjiK, 

£>2j+ 1-re-t  • 


^  ^  Amjk?(T  -  k)j  =  T2m+1,  as  f  =  T  £  M, 

fc=0  j 

From  ()1.25l)  we  can  obtain  the  following  identity  for  polynomials  in  T 
A  . _ \ _ t2j+1  I  2  A  (  — 


(1.26) 


j,  odd  k 


K,  J  2j  +  1  —  K 


B2j+i-KTK 


J  (2i  +  !)(?)' 

=  n2m+1,  asT->f. 

Taking  the  coefficient  of  n2-7-1-1  in  the  above  expression,  we  get  Amjin  =  (2m  +  1)(^),  and  taking 
the  coefficient  of  x2rf+1  for  an  integer  d  in  the  range  m/2  <  d  <  m  we  get  Am,d  =  0.  Taking  the 
coefficient  of  n2d+1  in  (11.261)  for  m/4  <  d  <  m/2  ,  we  get 


*-m,d 


(2<i+l)ff) 


■f  2(2  m  +  1) 


/  2m\  /  m 

\m  J  \2d  +  1 


i.e 


Am,d  —  (  1) 


m—  1 


(2m  +  1)! 


(-1  )m 
2m  —  2d 

1 


B2m—2d  —  0, 


~B2m—2d- 


d\d\m\(m  —  2d  —  1)!  m  —  d 
Continuing  similarly,  we  can  express  Arnj  for  each  integer  j  in  the  range  m/2s+1  <  j  <  m/2s 
(iterating  consecutively  s  =  1,2, ...)  via  the  previously  determined  values  of  d  <  j  as  follows 


Am, j  —  (2  j  +  1)  (  A  ^2  Am ) 


J 


,.d 


d  \ (-ly-1  D 
2j  +  1/  d-j  2d  2r 


d=2j+l 

The  same  formula  holds  also  for  m  =  0.  Note  that  the  m  in  above  sum  must  satisfy  m  >  2j  +  1  to 
return  a  nonzero  term  Amj. 

Definition  1.27.  We  define  here  a  generating  function  of  sequence  of  coefficients  Am  ]  as  follows 
0,  if  j  <  0  or  j  >  rn 

(2  j  +  1)  Cj)  Ed=2j+1  Am,d  {2f+1)  (~]-j  1  B2d—2j ,  if  0  <j<m 

.(2 j  +  l)(2/),  if  j  =  m 

Five  initial  rows  of  the  triangle  generated  by  Am  j,  j  >  0,  0  <  j  <  m  are 


4  — 

nm,j  • 
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m  =  0  1 

m  =  1  1 

m  =  2  1  0 

m  =  3  1  -14 

m  =  4  1  -120  0 

m  =  5  1  -1386  660 


Table  4.  Triangle  generated  by  Amj,  j  >  0, 


6 

30 

0  140 

0  630 

0  0  2772 


0  <  j  <  m,  see  definition  (1.27) 


Note  that  starting  from  row  m  >  11  the  terms  of  the  Itable  41  consist  of  fractional  numbers,  for 
example,  \  =  800361655623.6.  One  can  find  a  complete  list  of  the  numerators  and  denominators 
of  Am  j  in  OEIS  under  the  identifiers  A302971  and  A304042,  respectively,  see  |17].|18|.  Note  that 

=  22m+1  —  1. 

j 

As  we  have  found  a  recurrence  for  Amj  coefficients,  let’s  back  to  the  identity  CEO 

£  m  m 

E  E  Amjkj(T  -  ky  =  £(-i  )m-kum(e,  k )  •  Tk 

fc=l j= 0  k= 0 


to  find  a  recurrence  for  coefficients  Urn(i,  k)  for  every  m  >  1.  To  understand  the  nature  of  the 
coefficients  Um(T,  k),  express  the  polynomials  Um(£,  k)  in  terms  of  Amj  and  Bernoulli  numbers.  To 
do  so,  let’s  expand  the  binomial  (T  —  ky  in  the  l.h.s.  of  (11.11)  and  change  of  the  order  of  summation: 

(br'f-iy-v-' 


(1.28) 


k= 1  j= 0 


k= 1  j= 0  t= 0 

m  £  m  /  .\ 


=  EtTE  ) K/’i-i r‘. 


Now,  taking  the  coefficient  of  Tl  in  (I1.28P  gives: 


£  m 

fc=i  j=t 

From  this  formula  it  may  be  not  immediately  clear  why  Um(£,t)  represent  polynomials  in  T.  How¬ 
ever,  this  can  be  seen  if  we  change  the  summation  order  again  and  use  Faulhaber’s  formula  to 
obtain: 

(1.29)  Um(T,  t)  =  (-1)*"  E  (I)a^-  {Ziy  E  (23“t+1)B,T22-‘+1-'. 

i=t  ^  ^  + 1  ;=o  ^  ^ 

Introducing  k  =  2j  —  f  +  1  —  l  to  (11.291)  we  further  get  the  formula: 

um(T,t) . (-ir  E+1^E _r 

k — 1  j — t 

which  allows  easily  compute  the  coefficient  of  Tk  in  Um(T,t)  for  each  n.  In  above  formulae  we 
assume  that  B\  =  +i. 


Amj0-\-iy. 
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1.1.  Error  of  approximation.  Generally,  the  monomial  T2m+1,  m  €  N,  T  E  M  could  be 
approximated  by  as  follows 


lim 

T^ri 


i  m 

J2  Amtjkj{T  -  k)j 

=  lim 

T—ti 

m 

2(-l  )m-kUm(£,k)-Tk 

II 

II 

O 

,k= 0 

In  this  subsection  we  arrange  the  tables  with  values  of  corresponding  polynomials  Pm(£,T)  and 
monomials  T2m+1  to  rate  the  quality  of  the  approximation  of  T2m+1,  rn  >  0  by  (1.1).  We  begin 
from  the  case  m  =  1  and  T  =  3,  therefore,  the  approximating  polynomials  is  P\  3,  T  =  —81  +  36T. 
By  solving  the  equation  T3  —  36T  +  81  =  0  we  receive  the  following  roots  Ti  =  3,  T2  =  3.90833,  it 
follows  that  our  polynomial  approximates  the  monomial  T2m+1  in  neighborhood  of  T\  and  T2,  the 
following  table  contains  the  values  of  -Pi  (3,  T )  and  P3  in  neighborhood  of  the  root  T\  =  3 


T 

F(T)  =  T 3 

Pi(3,T)  =  -81  +  36T 

2.95 

25.67238 

25.2 

2.96 

25.93434 

25.56 

2.97 

26.19807 

25.92 

2.98 

26.46359 

26.28 

2.99 

26.73090 

26.64 

3.00 

27.00000 

27 

3.01 

27.27090 

27.36 

3.02 

27.54361 

27.72 

3.03 

27.81813 

28.08 

3.04 

28.09446 

28.44 

3.05 

28.37263 

28.8 

Table  5.  Table  of  values  of  F(T)  =  T3  and  Pi(3,T)  in  neighborhood  of  Ti  =  3. 


Below  we  show  graphically  the  intersection  of  of  F(T )  =  T3  and  Pi  (3,  T) 
F(T) 


Figure  4.  Intersection  of  F(T)  =  T3  and  Pl(3,T)  in  point  Ti  =  3. 
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As  it  is  mentioned  above  the  equation  T3  —  36 T  +  81  =  0  has  two  roots  T\  =  3,  T2  =  3.90833, 
let’s  discuss  the  case  T2  =  3.90833,  the  following  table  represents  the  values  of  P\  (3,  T)  and  T3  in 
neighborhood  of  the  root  T2  =  3.90833 


T 

F(T)  =  T6 

P\  (3,  T)  =  -81  +  36T 

3.86000 

57.51246 

57.96 

3.87000 

57.96060 

58.32 

3.88000 

58.41107 

58.68 

3.89000 

58.86387 

59.04 

3.90000 

59.31900 

59.4 

3.90833 

59.69991 

59.69988 

3.92000 

60.23629 

60.12 

3.93000 

60.69846 

60.48 

3.94000 

61.16298 

60.84 

3.95000 

61.62988 

61.2 

3.96000 

62.09914 

61.56 

Table  6.  Table  of  values  of  F(T)  =  T3  and  Pi(3,T)  in  neighborhood  of  =  3. 
Below  we  show  graphically  the  intersection  of  F(T)  =  T3  and  P\(3,T) 


F(T ) 


Next,  let’s  review  the  case  m  =  2  and  T  =  3,  the  approximating  polynomial  is  7-2(3,  T )  = 
2943  —  2160T  +  420T2.  By  solving  the  equation  T5  —  (2943  —  2160T  +  420T2)  =  0  we  receive  the 
following  roots  T\  =  3,  T2  =  3.40551,  T3  =  3.98704,  it  follows  that  our  polynomial  approximates 
the  monomial  T2m+1  in  neighborhood  of  T\,  T2  and  T3,  the  following  tables  contain  the  values  of 
7-2(3,  T )  and  T5  in  neighborhood  of  the  roots  T\  =  3,  T2  and  T3 
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T 

F(T)  =  T° 

P2( 3,  T)  =  2943  -  2160T  +  420T^ 

2.95 

223.41384 

226.05000 

2.96 

227.22628 

229.27200 

2.97 

231.09058 

232.57800 

2.98 

235.00728 

235.96800 

2.99 

238.97691 

239.44200 

3.00 

243.00000 

243.00000 

3.01 

247.07709 

246.64200 

3.02 

251.20872 

250.36800 

3.03 

255.39544 

254.17800 

3.04 

259.63780 

258.07200 

3.05 

263.93634 

262.05000 

Table  7.  Table  of  values  of  F(T)  =  T 5  and  P2(3,T)  in  neighborhood  of  T\  =  3. 


Below  we  show  graphically  the  intersection  of  F(T)  =  T'}  and  P-2(3,T) 
f(T) 


Next,  let’s  review  the  values  of  F(T)  =  T5  and  P2(3,T)  in  the  neighborhood  of  the  root  T2  = 
3.40551  as  follow  table  shows 


T 

F(T)  = 

P2{ 3,  T )  =  2943  -  2160T  +  420T^ 

3.35000 

421.91410 

420.45000 

3.36000 

428.24903 

427.03200 

3.37000 

434.65983 

433.69800 

3.38000 

441.14717 

440.44800 

3.39000 

447.71175 

447.28200 

3.40551 

458.04780 

458.04771 

3.41000 

461.07534 

461.20200 

3.42000 

467.87574 

468.28800 

3.43000 

474.75615 

475.45800 

3.44000 

481.71727 

482.71200 

3.45000 

488.75980 

490.05000 

Table  8.  Table  of  values  of  F(T)  =  T5  and  P2(3,T)  in  neighborhood  of  X2  =  3.40551. 
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Below  we  show  graphically  the  intersection  of  F(T )  =  T5  and  P2(3,T)  in  point  T2  =  3.40551 
f(T) 


Now,  let’s  review  the  values  of  F(T )  =  T'5  and  P2(3,T )  in  the  neighborhood  of  the  root  T2 
3.90833  as  follow  table  shows 


T 

F(T)  = 

P2{ 3,  T)  =  2943  -  2160T  +  420T^ 

3.93000 

937.48160 

941.05800 

3.94000 

949.46970 

952.51200 

3.95000 

961.58012 

964.05000 

3.96000 

973.81381 

975.67200 

3.97000 

986.17170 

987.37800 

3.98000 

998.65472 

999.16800 

3.98704 

1007.51835 

1007.51854 

4.00000 

1024.00000 

1023.00000 

4.01000 

1036.86416 

1035.04200 

4.02000 

1049.85728 

1047.16800 

4.03000 

1062.98034 

1059.37800 

Table  9.  Table  of  values  of  F(T)  =  T5  and  P2(3,T)  in  neighborhood  of  X3  =  3.98704. 


Below  we  show  graphically  the  intersection  of  F(T)  =  T 5  and  p2(3,T)  in  point  T3  =  3.98704 
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f(T) 


Now,  let’s  keep  our  attention  to  the  case  m  =  3  and  1  =  2,  the  approximating  polynomial  is 
P3( 2,  T)  =  -9028  +  13818T  -  7140T2  +  1260T3.  By  solving  the  equation  T7  -  (-9028  +  13818T  - 
7140T2  +  1260T3)  =  0  we  receive  the  following  roots  Ti  =  2,  T2  =  2.8605,  T3  =  2.98955,  it  follows 
that  our  polynomial  approximates  the  monomial  T‘  in  neighborhood  of  Tf,  T2  and  T3,  the  following 
tables  contain  the  values  of  P2(3,T)  and  T5  in  neighborhood  of  the  roots  T\,  T2  and  T3,  note  that 
we  have  arranged  the  roots  T2  and  T3  in  single  table 


T 

F(T)  =  T‘ 

7-3(2,  T)  =  -9028  +  13818T  -  7140T2  +  1260T3 

1.96 

111.12007 

113.47136 

1.97 

115.14990 

117.00398 

1.98 

119.30436 

120.59792 

1.99 

123.58664 

124.26074 

2.00 

128.00000 

128 

2.01 

132.54776 

131.82326 

2.02 

137.23333 

135.73808 

2.03 

142.06015 

139.75202 

2.04 

147.03177 

143.87264 

2.05 

152.15178 

148.1075 

2.06 

157.42385 

152.46416 

Table  10.  Table  of  values  of  F(T)  =  T7  and  P3(2,T)  =  -9028  +  13818T-7140T2  +  1260T3 
in  neighborhood  of  Ti  =  2. 


Below  we  arrange  a  plot  of  intersection  ofP(T)  =  T‘  andP3(2,T)  =  —  9028  +  13818T  —  7140T2  + 
1260T3  in  Tx  =  2 
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f(T) 


Let’s  show  the  values  of  the  F(T)  =  T~  and  polynomial  P3(2,  T )  =  —9028  +  13818T  —  7140T2 
1260T3  in  neighborhood  of  two  other  roots,  T\  =  2.86046,  T2  =  2.98955 


T 

F(T)  =  T‘ 

P3( 2,  T)  =  -9028  +  13818T  -  7140T2  +  1260T3 

2.81000 

1383.38875 

1379.35766 

2.82000 

1418.22050 

1415.09168 

2.83000 

1453.80129 

1451.52962 

2.84000 

1490.14449 

1488.67904 

2.85000 

1527.26367 

1526.5475 

2.86046 

1566.94022 

1566.940217 

2.87000 

1603.88519 

1604.47178 

2.88000 

1643.41563 

1644.54272 

2.89000 

1683.77827 

1685.36294 

2.90000 

1724.98763 

1726.94 

2.91000 

1767.05848 

1769.28146 

2.92000 

1810.00577 

1812.39488 

2.93000 

1853.84467 

1856.28782 

2.94000 

1898.59054 

1900.96784 

2.95000 

1944.25897 

1946.4425 

2.96000 

1990.86576 

1992.71936 

2.97000 

2038.42692 

2039.80598 

2.98000 

2086.95867 

2087.70992 

2.98955 

2134.25217 

2134.252175 

3.00000 

2187.00000 

2186 

3.01000 

2238.54314 

2236.40126 

3.02000 

2291.12403 

2287.65008 

3.03000 

2344.76002 

2339.75402 

3.04000 

2399.46868 

2392.72064 

Table  11.  Table  of  values  of  F(T)  =  T7  and  P3(2,T)  =  -9028  +  13818T-7140T2  +  1260T3 
in  neighborhood  of  Ti  =  2.8605,  T3  =  2.98955. 


Below  we  show  graphically  the  intersection  of  F(T )  =  T‘  and  P3(2,T)  =  —9028  +  13818T 
7140T2  +  1260T3  in  point  T3  =  3.98704  in  points  T2  =  2.8605,  T3  =  2.98955 
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3.  Conclusion 

In  this  manuscript  we  have  derived  and  discussed  the  polynomials  Pm(l,  T),  such  that  converge  to 
the  monomial  F2m+1,  m  E  No  as  T  — >  £,  £  €  N.  Polynomial  Pm(£,T)  gives  weak  approximation  of 
the  monomial  T2m+1  mainly  in  neighborhood  of  the  solution  of  the  equation  T2m+1  —  Pm(£,  T )  =  0. 

4.  Supplementary  Files 

We  provide  the  following  supplementary  files  to  support  our  study: 

•  https:/ /kolosovpetro.github.io/arxiv_1603_02468/identity_l_l.txt  -  Mathematica  program,  im¬ 
plementation  the  l.h.s  of  identity  (11.11).  that  is  T2m+1  =  Ylk=i  Xq=o  Am,j^(T  —  k)\  £  = 
TeN. 

•  https:/ /kolosovpetro.github.io/arxiv_1603_02468/identity_l_l_r_h_s.txt  -  Mathematica  program, 

implementation  the  r.h.s  of  identity  (11.11).  that  is  T2m+1  =  l)m_fcfym(£,  k)  ■  Tk ,  l  = 

TeN. 

•  https://kolosovpetro.github.io/arxiv_1603_02468/u_coefficients_in_row.txt  -  Mathematica  pro¬ 
gram,  lists  the  coefficients  Um(£,  k)  in  (11,11)  for  given  m.  (By  default  rn  =  2) 
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